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Abstract 

In this note we present an extension of our previous analysis of some pre- 
liminary experimental results of the DAMA/Nal Collaboration which might 
be indicative of a yearly modulation effect. Here we present a direct way for 
obtaining from the experimental data the relevant cosmological implications 
for relic neutralinos. We find that some of the configurations singled out by 
the DAMA/Nal results would have cosmological properties compatible with 
a neutralino as a dominant component of cold dark matter (on the average in 
the Universe and in our galactic halo). 
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Some recent results of the DAMA/Nal Collaboration on direct WIMP detection seem 
to indicate that a yearly modulation effect might be present. As stressed in Ref. [|l], a much 
larger statistics is necessary, before any claim of a real effect can be made; and, actually, a 
further collection of data is under way at the large-mass, low-background Nal(Tl) detector 
at the Gran Sasso Laboratory [I]]. Meanwhile, the preliminary results are so intriguing 
to prompt an analysis of their possible implications for particle physics and cosmological 
properties. 

In Ref. 0, we have analyzed some possible consequences for neutralino properties implied 
by the results of Ref. under the hypothesis that these are indicative of a yearly modulation 
effect due to relic neutralinos. We have shown that a number of neutralino configurations, 
compatible with all current experimental constraints, might produce a modulation effect at 
the level of the one possibly detected by the DAMA/Nal Group. Some appealing features of 
these configurations (or at least of a subset of these) are: i) sizeable cosmological density Q, 
ii) detectability of indirect signals at neutrino telescopes, iii) detectability of some related 
properties at accelerators. 

In our previous paper [Q] the comparison of theoretical evaluations with experimental 
data was performed in terms of a plot in the plane £&icailr° n ^ vs - m x> where m x is the 
neutralino mass, o"g°^ 1 r eon ' ) is the scalar elastic cross section off nucleon and £ = p x j pi is the 
fractional amount of local neutralino dark matter density p x with respect to the total local 
dark matter density p\. The results are reported in Fig. 1 (which is taken from Ref. 0). The 
closed contour line denotes the region (R m ) singled out at 90% C.L. by the data of Ref. M. 
The open curve denotes the 90% C.L. upper bound of C,cr^i^ on \ as obtained from the total 
counting rates of Ref. ||. In extracting both the open and the closed contour lines from the 
experimental data, the values of some astrophysical parameters relevant for the event rates 
at the detector have to be chosen. These parameters are: the root mean square velocity v rms 
of the neutralino Maxwellian velocity distribution in the halo, the neutralino escape velocity 
v esc in the halo, the velocity v & of the Sun around the galactic centre. The values adopted 
in Fig. 1 refer to the median values of these parameters in their experimentally allowed 
ranges M, namely: t> rms = 270 Km s _1 , v esc = 650 Km s -1 , v & = 232 Km s -1 . We remark 
that the values of £&icailr° n \ extracted from experimental data, are not very sensitive to the 
specific values adopted for the various velocities within their physical ranges. The values of 
£ a scaiar !0n ' ) plotted in Fig. 1 are normalized to the value p\ = 0.5 GeV cm -3 . 

Together with the experimental results, Fig. 1 shows our evaluations of the quantity 
£ CT scaiax° inside the Minimal Supersymmetric Standard Model. For a brief discussion of the 
model we are employing, the reader is addressed to Ref. 0. Here we only recall that the 
supersymmetric parameter space, considered here and in Ref. ||, consists of six independent 
parameters: M 2 , p, tan/3, m^, m , A. We vary these parameters in the following ranges: 
10 GeV < M 2 < 500 GeV (21 steps over a linear grid); 10 GeV < \p\ < 500 GeV (21 
steps, linear grid); 65 GeV < rriA < 500 GeV (15 steps, logarithmic grid); 100 GeV < m < 
500 GeV (5 steps, linear grid); — 3 < A < +3 (5 steps, linear grid); 1.01 < tan/? < 50 (15 
steps, logarithmic grid). 

The supersymmetric parameter space is constrained by the experimental limits obtained 
from accelerators on supersymmetric and Higgs searches and from the b — > s + 7 radiative 
decay. For a complete set of references, see Ref. [[|. In addition to the experimental limits, 
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the parameter space is further constrained by the requirement that the neutralino is the 
lightest super symmetric particle, since the neutralino is the dark matter candidate under 
investigation in the present analysis. Finally, the regions of the parameter space where the 
neutralino relic abundance exceeds the cosmological bound, i.e. Q x h 2 > 1, are also excluded. 

Our analysis in Ref. was performed employing, for the evaluation of £, the usual 
rescaling recipe || which compares the amount of neutralino relic density Q x h 2 to the mini- 
mal value for the total amount of dark matter in the Universe compatible with observational 
data and with large-scale structure calculations. The recipe works as follows: when Q x h 2 
is larger than (f2/i 2 ) min we simply put £ = 1; when Q x h 2 turns out to be less than (Qh 2 ) min , 
and then the neutralino may only provide a fractional contribution Q x h 2 /(Qh 2 ) min to Qh 2 , 
we take £ = Q x h 2 /(Qh 2 ) m i n . 

Besides the empirical character of the rescaling procedure, it has to be noted that, as was 
stressed in Ref. @, the value of (f2/i 2 ) min is rather arbitrary, possibly in the range 0.03 < 
(Qh 2 ) min < 0.3. To be definite, in Ref. we reported our results for the representative 
value (flh 2 ) min = 0.03, and this is the value used in the scatter plot of Fig. 1. 

In this note we present an extension of our previous investigation, which is meant to ob- 
tain from the experimental data the relevant cosmological implications for relic neutralinos, 
in the most direct way without any use of rescaling for the neutralino local density. 

The procedure we adopt is the following: 1) We evaluate cr s ^ 1 ^ 1 r eon ' ) and Vt x h 2 by vary- 
ing the supersymmetric parameters over the previously defined grid. 2) For any value of 
[Px CT scaiar ]expt = Pi [C^scajar 1 expt compatible with the experimental region R m we calculate 
p x as given by p x = Ipx^e^^ 11 ]expt/^^^° an d restrict the values of m x to stay inside the 
region R m . 3) The results are displayed in a scatter plot in the plane p x vs. Q x h 2 . 

Examples of our results are given in Figs. 2-4 for a few experimentally allowed values 
of £cr^ 1 ^ 1 r eon ' ) , taken from Fig. 1 (remember that the normalization value p\ = 0.5 GeV cm -3 
is used in this figure). Figs. 2,3,4 refer to [C^scai^T^lexpt = 6 • 10" 9 nbarn, 2 • 10~ 9 nbarn 
and 1 ■ 10~ 9 nbarn, respectively. The two horizontal lines denote the physical range 0.1 
GeV cm" 3 < pi < 0.7 GeV cm" 3 for the local density of non-baryonic dark matter. This 
(rather generous) range has been established by taking into account possible effects on pi 
due to a flattened dark halo || and a recent derivation of the contribution of baryons to 
the local dark matter from microlensing data. The solid vertical lines denote the physical 
band for Qh 2 : 0.03 < Qh 2 < 1, and the two dashed lines give the favored band for the cold 
dark matter contribution to fl: 0.1 < (H/t 2 )cDM < 0.3 f|. The two tilted dot-dashed lines 
denote the band where linear rescaling procedure for the local density is usually applied. 
The upper dot-dashed line refers to a rescaling with (Qh 2 ) min = 0.03, the lower one to the 
value (Qh 2 ) min = 0.3. With the aid of this kind of plot we can classify the supersymmetric 
configurations belonging to region R m into various categories. Configurations whose repre- 
sentative points fall above the maximum value p x = 0.7 GeV cm -3 have to be excluded (we 
remind that those providing an fl x h 2 > 1 are already disregarded from the very beginning). 
Among the allowed configurations, those falling in the region inside both the horizontal and 
solid vertical lines (this region is called A hereafter) are very appealing, since they would 
represent situations where the neutralino could have the role of a dominant cold dark matter 
component; even more so, if the representative points fall in the subregion (B) inside the 
vertical band delimited by dashed lines. Configurations which fall inside the band delimited 
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by the tilted dot-dashed lines denotes situations where the neutralino can only provide a 
fraction of the cold dark matter both at the level of local density and at the level of the 
average Q. Configurations above the upper dot-dashed line and below the upper solid hori- 
zontal line would imply an unlikely special clustering of neutralinos in our halo as compared 
to their average distribution in the Universe. 

It is worth noticing a few important properties of the scatter plots shown in Figs. 2-4: 

1) The scatter plots display a correlation between p x and Q x h 2 . This feature is expected 
on the basis of the following properties: i) Q x h 2 is roughly inversely proportional to 
the neutralino pair annihilation cross section, ii) at fixed [C^icaiaT lexpt, Px * s inversely 
proportional to , iii) the annihilation cross section and cr^^ on ^ are usually 
correlated functions (i.e., they are both increasing or decreasing functions of the su- 
persymmetric parameters). 

2) The regions A and B become more populated in configurations as the experimental 
value [^Cs^iaT" 1 lexpt decreases (i.e., if one proceeds from the case shown in Fig. 2 to the 
one in Fig. 4). This feature follows from the fact that cr^^ on ^ is bounded from above 
by accelerator limits (mainly because of lower bounds on Higgs masses); this implies 
for p x a lower bound, which however is less stringent at lower values of [^scalar ] e xpt- 

3) It is remarkable that, as a consequence of point 2), some neutralino configurations 
populate the most interesting region, i.e. region B. This occurs for [£crg™k!r ]expt < 
6 ■ 10~ 9 nbarn. 

In what follows we concentrate on configurations falling in region A or in the subregion B. 
For definiteness, we choose the representative value [C^scaiaT^lexpt = 2-1CT 9 nbarn. First, we 
check whether some of these neutralino compositions are excluded by indirect detection of 
WIMPs, by measurements of up-going muon fluxes due to the neutrino outburst produced by 
neutralino pair annihilation in the Earth or in the Sun. In Fig. 5 we show how our calculated 
fluxes from the Earth $^ arth compare with the current experimental upper bounds |),[I0 . 
We notice that some configurations are actually excluded by the present limits on $^ arth . 
From our calculations, it turns out that the present upper bounds on the up-going muon 
fluxes from the Sun constrain the analyzed configurations only marginally. 

Finally, we provide in Figs. 6,7 some information about important parameters referring 
to the configurations of region A and of subregion B. These figures also show what may be 
the exploration potential for these configurations at LEP2 and Tevatron. 

In conclusion, we remark that our present analysis shows, in a way independent of the 
rescaling procedure on local DM density, that some of the configurations singled out by the 
DAMA/Nal results would have cosmological properties compatible with a neutralino as a 
dominant component of cold dark matter (on the average in the Universe and in our galactic 
halo) . 
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FIGURE CAPTIONS 



Figure 1 - The scalar neutralino-nucleon cross section (T^iar ■> multiplied by the 
factor £, is plotted versus the neutralino mass m x . The closed contour delimits the 
region singled out at 90% C.L. by the data of Ref. fflj. The open curve denotes 
the 90% C.L. upper bound obtained from the total counting rates of Ref. ||. The 
scatter plot represents the theoretical prediction for Create > calculated within the 
MSSM scheme. Only configurations with /i > are displayed. The cosmological and 
astrophysical parameters are set to the following values: v rms = 270 Km s" 1 , t> esc = 650 
Km s"\ v Q = 232 Km s~\ p, = 0.5 GeV cm" 3 and (nh 2 ) min = 0.03. 

Figure 2 - The neutralino local density p x , calculated for [^scaiaT" 1 ]expt = 6 ■ 10~ 9 
nbarn, is plotted versus the neutralino relic abundance Q x h 2 . For the value of 

[C^scaiaf 3 ]expt employed here, the neutralino mass is restricted to the range 30 GeV 
< m x < 155 GeV, as obtained from the closed contour in Fig. 1. The two horizontal 
lines denote the physical range for the local density of non-baryonic dark matter. The 
two solid vertical lines denote the physical band for Qh 2 , and the two dashed lines 
give the preferred band for the cold dark matter contribution to Q. The two tilted 
dot-dashed lines denote the band where linear rescaling procedure for the local density 
is usually applied. The region inside both the two horizontal lines and the two vertical 
solid (dashed) lines is defined as region A (region B). 

Figure 3 - The neutralino local density p x , calculated for [^o"^^ 3 ] eX pt — 2 ■ 10~ 9 
nbarn, is plotted versus the neutralino relic abundance Q x h 2 . The neutralino mass is 
restricted to the range 30 GeV < m x < 115 GeV, as obtained from the closed contour 
in Fig. 1. The solid, dashed and dot-dashed lines, as well as region A and region B, 
are defined as in Fig. 2. 

Figure 4 - The neutralino local density p x , calculated for [£cscaia^° ] eX pt = 1 • 10~ 9 
nbarn, is plotted versus the neutralino relic abundance Q x h 2 . The neutralino mass is 
restricted to the range 40 GeV < m v < 85 GeV, as obtained from the closed contour 
in Fig. 1. The solid, dashed and dot-dashed lines, as well as region A and region B, 
are defined as in Fig. 2. 

Figure 5 - Flux of up-going muons $^ arth as a function of the neutralino mass m x , 
calculated for the neutralino configurations belonging to region A (dots) and to region 
B (crosses) of Fig. 3 (therefore referring to [^icX^^cxpt = 2 ■ 10~ 9 nbarn and 40 
GeV < m x < 115 GeV). The solid and dashed lines represent the experimental 90% 
C.L. upper bounds of Ref. || and Ref. fl(|, respectively. 



Figure 6 - The configurations compatible with region A (dots) and region B (crosses) 
of Fig. 3 are plotted in the m x -tan (3 plane. Configurations which provide a flux of 
up-going muons exceeding the experimental upper limits of Ref. flJToH are dropped. 



The dark region on the left side is excluded by current LEP data |[T|. The region on 



the left of the vertical solid line will be accessible to LEP at y/s = 192 GeV ]T2]. The 



region on the left of the vertical dashed line will be explorable at TeV33 [O 
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Figure 7 - The configurations compatible with region A (dots) and region B (crosses) 
of Fig. 3 are plotted in the m^-tan (3 plane. Configurations which provide a flux of up- 
going muons exceeding the experimental upper limits of Ref. P,|T0| are dropped. The 
dark regions are excluded by current LEP searches or by theoretical arguments. 
The region on the left of the solid line will be accessible to LEP at y/s = 192 GeV, 
with a luminosity of 150 pb _1 per experiment Il2l . 
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